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Abstract. The influence of the gravity acceleration on the regularized energy- momentum 
tensor of the quantized electromagnetic field between two plane parallel conducting plates 
is derived. We use Fermi coordinates and work to first order in the constant acceleration 
parameter. A new simple formula for the trace anomaly is found to first order in the 
constant acceleration, and a more systematic derivation is therefore obtained of the theo- 
retical prediction according to which the Casimir device in a weak gravitational field will 
experience a tiny push in the upwards direction. 
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1. Introduction 

An important property of quantum electrodynamics is that suitable differences of zero- 
point energies of the quantized electromagnetic field can be made finite and produce mea- 
surable effects such as the tiny attractive force among perfectly conducting parallel plates 
known as the Casimir effect [1]. This is a remarkable quantum mechanical effect that 
makes itself manifest on a macroscopic scale. For perfect reflectors and metals the Casimir 
force can be attractive or repulsive, depending on the geometry of the cavity, whereas for 
dielectrics in the weak-reflector approximation it is always attractive, independently of the 
geometry [2] . The Casimir effect can be studied within the framework of boundary effects 
in quantum field theory, combined with zeta-function regularization or Green-function 
methods, or in more physical terms, i.e. on considering van der Waals forces [3] or scatter- 
ing problems [4] . Casimir energies are also relevant in the attempt of building a quantum 
theory of gravity and of the universe [5]. 

For these reasons, in Ref. [6] we evaluated the force produced by a weak gravitational 
field on a rigid Casimir cavity. Interestingly, the resulting force was found to have opposite 
direction with respect to the gravitational acceleration; moreover, we found that the current 
experimental sensitivity of small force macroscopic detectors would make it possible, at 
least in principle, to measure such an effect [6]. More precisely, the gravitational force 
on the Casimir cavity might be measured provided one were able to use rigid cavities 
and find an efficient force modulation method [6]. Rigid cavities, composed by metal 
layers separated by a dielectric layer, make it possible to reach separations as small as 
5-j-10 nm and allow to build multi-cavity structures, made by a sequence of such alternate 
layers. If an efficient modulation method could be found, it would be possible to achieve 
a modulated force of order 10~ 14 N in the earth's gravitational field. The measure of such 
a force, already possible with current small-force detectors on macroscopic bodies, might 
open the way to the first test of the gravitational field influence on vacuum energy [6] . In 
Ref. [6], calculations were based on simple assumptions and the result can be viewed as a 
reasonable " 'first order 1 '' generalization of T Ml/ from Minkowski to curved space-time. The 
present paper is devoted to a deeper understanding and to more systematic calculations 
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of the interaction of a weak gravitational field with a Casimir cavity. To first order in our 
approximation the former value of the force exerted by the field on the cavity is recovered. 
But here we also find a trace anomaly for the energy-momentum tensor. 

We consider a plane-parallel Casimir cavity, made of ideal metallic plates, at rest in 
the gravitational field of the earth, with its plates lying in a horizontal plane. We evaluate 
the influence of the gravity acceleration g on the Casimir cavity but neglect any variation 
of the gravity acceleration across the cavity, and therefore we do not consider the influence 
of tidal forces. The separation a between the plates is taken to be much smaller than the 
extension of the plates, so that edge effects can be neglected. We obtain a perturbative 
expansion of the energy-momentum tensor of the electromagnetic field inside the cavity, 
in terms of the small parameter e = 2ga/c 2 , to first order in e. For this purpose, we use a 
Fermi [7,8] coordinates system (t,x,y,z) rigidly connected to the cavity. The construction 
of these coordinates involves only invariant quantities such as the observer's proper time, 
geodesic distances from the world-line, and components of tensors with respect to a tetrad 
[8]. This feature makes it possible to obtain a clear identification of the various terms 
occurring in the metric. In our analysis we adopt the covariant point-splitting procedure 
[9,10] to compute the perturbative expansion of the relevant Green functions. Gauge 
invariance plays a crucial role and we check it up to first order by means of the Ward 
identity. We also evaluate the Casimir energy and pressure, and in this way we obtain a 
sound derivation of the result in Ref. [6] , according to which the Casimir device in a weak 
gravitational field will experience a tiny push in the upwards direction. 

Use is here made of mixed boundary conditions on the potential plus Dirichlet con- 
ditions on ghost fields. With our notation, the z-axis coincides with the vertical upwards 
direction, while the (x, y) coordinates span the plates, whose equations are z = and 
z = a, respectively. The resulting line element for a non-rotating system is therefore [7] 

ds 2 = -c 2 (l + e-) dt 2 + dx 2 + dy 2 + dz 2 + 0(\x\ 2 ) = r)^dx*dx v - e-c 2 dt 2 , (1.1) 

where rj^ is the flat Minkowski metric diag(— 1, 1, 1, 1). 



3 



2. Green Functions 



For any field theory, once that the invertible differential operator Uij in the functional 
integral is given, the corresponding Green functions satisfy the condition (we use hereafter 
the DeWitt condensed-index notation) 

Ua(P k = -*<* (2.1) 

and are boundary values of holomorphic functions. The choice of boundary conditions will 
determine whether we deal with advanced Green functions G +jfc , for which the integration 
contour passes below the poles of the integrand on the real axis, or retarded Green functions 
for which the contour passes instead above all poles on the real axis, or yet other 
types of Green functions. Among these, a key role is played by the Feynman Green function 
G^p, obtained by choosing a contour that passes below the poles of the integrand that lie 
on the negative real axis and above the poles on the positive real axis. If one further 
defines the Green function [11] 

G jk = l(G+i k + G~ jk ), (2.2) 

one finds in stationary backgrounds (for which the metric is independent of the time 
coordinate, so that there exists a timelike Killing vector field) that the Feynman Green 

function has a real part equal to G^ , and an imaginary part equal to the Hadamard 
function iP fc , i.e. 

W k (x, x') = -2i \G J F h (x, x') - G jk (x, x')] ■ (2.3) 

This relation can be retained as a definition when the background is nonstationary; in such 
a case, however, H^ k (x,x') is generally no longer real [11]. 

In particular, the photon Green function G\ u > in a curved spacetime with metric g^ u 
solves the equation [12] 

yf=gPf(x)G Xv .= giiV 8(x,x'). (2.4) 
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The wave operator P x results from the gauge-fixed action with Lorenz [13] gauge-fixing 
functional <&l{A) = VM M , and having set to 1 the gauge parameter of the general theory, 
so that 

P lt X (x) = -8tUx + Ri(x), (2.5) 

where [] x = g a/3 V ' a V ' p{x) . Since we need the action of the gauge-field operator P^{x) 
on the photon Green functions, it is worth noticing that 



(2.6) 



The Christoffel coefficients for our metric (1.1) read 



r^ 7 = \g a5 (gsf>„ + gw - gw) = -\\ {v a °S°^ + V a °S% - rf<*8%) . (2.8) 

Since the connection coefficients, to first order in e, are constant, we realize that the 
Ricci curvature tensor vanishes to this order. On expanding (this is, in general, only an 
asymptotic expansion) 

Gv~G2+eGW+0(e 2 ), (2.9) 



we get 



so that 



QafiXv' = dadpGxu' — d 



r t* Mo) 

L fiX^/iv' 



-K^G^-r^G^,, 



and eventually 



r (o) , r (i) 



(i) _ r v r»(°) 



~ L aP^Xv',11 L aX (j fiu',/3 



e-SS6gd a dpG^. 
ct 



(2.10) 



(2.11) 



(2.12) 
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We therefore get, to first order in e, 



\-\o r (o) _ AO) 



(2.13) 



\—\o r (i) - 



(1) 



(2.14) 



where 



J { ^ = -r}^8{x,x'), (2.15) 
eJ -' " I 6 (¥ + ^) ^ + 2 ^ r ^,p + »T r^G$ , T " ^Coo, (2-16) 

with n° = ^dadp = -dl + d 2 x + d 2 y + d 2 z . 

To fix the boundary conditions we note that, on denoting by E t and H n the tangen- 
tial and normal components of the electric and magnetic fields, respectively, a sufficient 
condition to obtain 



0, H„ 



0, 



(2.17) 



on the boundary S of the device, is to impose Dirichlet boundary conditions on 

A (x),A 1 (x),A 2 (x) 

[14] at the boundary z = 0, z = a. The boundary condition on As is determined by 
requiring that the gauge-fixing functional, here chosen to be of the Lorenz type, should 
vanish on the boundary. This implies 



A^\ s = =► A%\ s = (g™d 3 A 3 - gn%,A*)\ 8 = • 



(2.18) 



To first order in e, these conditions imply the following equations for Green functions: 



0, li = 0,1,2, W, 



= 0, Vi/, 



= 0, /i = 0,1,2, W, 



(2.19) 

(2.20) 
(2.21) 
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• ^.(0) 



, Vi/, (2.22) 

s 



s 2a 3i/ ' 

hence we find that the third component of the potential satisfies homogeneous Neumann 
boundary conditions to zeroth order in e and inhomogeneous boundary conditions to first 
order. 

Since jj^J, is diagonal, by virtue of the homogeneity of the boundary conditions, G^, 
turns out to be diagonal. On the contrary, jj^J, has two off-diagonal contributions: Jq\' 
and J$q\ so that G^, is non-diagonal. Let us write down explicitly the expressions for 
the various components of J^J,, i.e. 

a ^oJ = 2&(x,x ) — zGqq, Q0 + 2^oo^',3' (2.23) 

a ^03' = _ ^33',0> (2.24) 

ajff = Z -S(x, x>) - zG?l 00 - \Gf^ 3 , (2.25) 

a ^22' = 2^( X ' X ) ~~ Z< ^22', 00 _ 2^22', 3' (2.26) 
a ^33' = 2~d(x,X ) — zG^ / QQ — 2 ^33', 3' (2.27) 

a ^io' = _ ^oo',o- (2.28) 

Now we are in a position to evaluate, at least formally (see below), the solutions to zeroth 
and first order, and we get 

r<(o) „ f ducPk iul (t-t')+ik_ L -(x_ L -x'j ) n /„ „/\ (o 0Q \ 

G \u' = V\v' j e ± '9d,n{z,z), (2.29) 

having defined 

, . sin /c(az<) sin n(a — z>) , 

g D (z,z'-K) = — , 0<z,z'<a, 2.30 

k sin KCl 

, . cos «(az<) cos K,(a — z>) , 

g N (z, z'; k) = — , < z, z' < a, 2.31 

Ksmna 
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where D, N stand for homogeneous Dirichlet or Neumann boundary conditions, respec- 
tively, z > (z<) are the larger (smaller) between z and z', while k± has components (k x , k y ), 
x± has components (x,y), k = \Ju) 2 — /c 2 , and 

C = / ^e--^-*')^-^-^)^, (2.32) 

where the $ components different from zero are written in Ref. [15]. A scalar field satisfies 
the same equations of the 22 component of the gauge field, hence we do not write it 
explicitly. In the following we will write simply G^ v > and G for the Green function of the 
gauge and ghost field, respectively. 

We should stress at this stage that, in general, the integrals defining the Green func- 
tions are divergent. They are well defined until x ^ x', hence we will perform all our 
calculations maintaining the points separated and only in the very end shall we take the 
coincidence limit [16]. We have decided to write the divergent terms explicitly 

so as to bear them in mind and remove them only in the final calculations by hand, instead 
of making the subtraction at an earlier stage. 

Incidentally, we note that these Green functions satisfy the Ward identity 

G V,n + G -y = °> G V' + G ^ = °> ( 2 - 33 ) 

to first order in e so that, to this order, gauge invariance is explicitly preserved (the check 
being simple but non-trivial). 

3. Energy-Momentum Tensor 

By virtue of the formulae of Sec. 2 we get, from the asymptotic expansion T^ v > ~ Tjfj, + 

fr<:;+o(a, 
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where (h is the Hurwitz ^-function £u(x,/3) = Yl^Lo( n + 13) x . On taking the limit 
z' — > z + we get 

}™J^') = + < lim -3^) diag(-l, 1, 1, -3), (3.2) 

where the divergent term as z' — > z can be removed by subtracting the contribution of 
infinite space without bounding surfaces [1], and in our analysis we therefore discard it 
hereafter. The renormalization of the energy-momentum tensor in curved spacetime is 
usually performed by subtracting the {T^ u ) constructed with an Hadamard or Schwinger- 
DeWitt two-point function up to the fourth adiabatic order [9,17]. In our problem, however, 
as we work to first order in e, we are neglecting tidal forces and therefore the geometry of 
spacetime in between the plates is flat. Thus, we need only subtract the contribution to 
the energy momentum tensor that is independent of a, which is the standard subtraction 
in the context of the Casimir effect in flat spacetime. 
In the same way we get, to first order in e: 

lim (T( 1 ^') = diag(T( 1 )° ,TW 11 ,T( 1 ) 22 ,T( 1 ) 33 )+ lim diagf - z'/n 2 (z - z')\ 0, 0, o), 

(3.3) 



where 



T (i)oo = ^ , , ^cot(^)csc 2 (^) 

1200 a 3 3600 a 4 240 a 3 ' 1 ' ' 

T q)ii 7 rcot(^)csc 2 (^) 

3600 a 3 1800 a 4 120 a 3 ' 1 ' ' 

T (i)22 = T (i)n (3 _ 6) 

(3.7) 



T (i)33 = (?r 2 (a-2z)) 
720 a 4 



By virtue of the Ward identities for quantum electrodynamics, here checked up to first 
order in e, the gauge-breaking part of the energy-momentum tensor is found to be minus 
the ghost part, hence we compute the second only. 
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4. Casimir Energy and Force 

To compute the Casimir energy we must project the energy-momentum tensor along the 
unit timelike vector u with covariant components = (y/— goo, 0, 0, 0) to obtain p = 
(Ti* v )upU v , so that 



(^) 



TT 2 e / TT 2 : n 2 z 7rcot(^)csc 2 (^) 



720a 4 a V 1200 a 3 3600 a 4 240 a 3 



7T 



2 (7 / 7T 2 n 2 z n cot(^) csc 2 (^, . 



720a 4 c 2 V 1200 a 3 900 a 4 240 a 3 

where in the second line we have substituted e by its expression in terms of g. Thus, the 
energy stored in the Casimir device is found to be 

where A is the area of the plates, d 3 E is the three-volume element of an observer with 
four-velocity u M , Eq. (4.2) is expressed through a principal-value integral, and we have 
reintroduced h and c. 

In the same way, the pressure on the plates is given by 

^/ i-\ \ 71-2 / 2 qa\ „. . 7T 2 ac / 2 qa\ . , „. 

P(S = °» = 240^ I 1 + S?) ' P( * = = - 240^ I 1 - 8?) ' (4 ' 3) 

so that a net force pointing upwards along the z-axis is obtained, in full agreement with 
Eq. (8) in Ref. [6], with magnitude 

^ 7T 2 Ahq , jt . 

F = f . (4.4) 

180 ca 3 V ; 

The reader may wonder whether the pressure on the outer faces of the cavity may alter this 
result. A simple way to answer this question is to imagine that our cavity is included into 
a surrounding cavity on both sides. On denoting by b the common separation between 
either plates of the original cavity and the nearest plate of the surrounding cavity, and 
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assuming that b is such that a/b « 1, but still small enough so as to obtain gb/c 2 « 1, 
we see from Sec. 3 that the outer pressure on both plates of the original cavity includes the 
same divergent contribution which acts from within plus a finite contribution that becomes 
negligible for a/b « 1. To sum up, the divergent contributions to the pressure from the 
inside and the outside of either plate cancel each other exactly, and one is left just with 
the finite contribution from the inside, as given in Eq. (4.3). 

As a check of the result, it can be verified that the energy- momentum tensor is covari- 
antly conserved to first order in e. To this order, the covariant conservation law implies 
the following conditions: 

e° : (Tt )""),, = 0, (4.5) 

e : (T^)j = (i = 0, 1, 2), i ((T^ 00 ) + (T^ 33 )) + (T (1)33 ), 3 = 0, (4.6) 

that are indeed satisfied. Moreover, from the previous expressions of the energy-momentum 
tensor the following trace anomaly r is obtained: 

The volume integral of this density exists as a principal-value integral and is given by 

fr^=^%A. (4.8) 
J 360 car 

The global, integrated form (4.8) of the trace anomaly Eq. (4.7) is the new result with 
respect to the analysis in Ref . [6] . It tends to zero at large separation a between the plates. 
This trace anomaly is therefore caused by the presence of the boundaries, and then is of a 
different nature from the usual trace anomaly which is encountered in curved spacetimes 
without boundaries, which depends on the Riemann curvature [11,16]. 

5. Concluding Remarks 

To the best of our knowledge, the analysis presented in this paper represents the first 
study of the energy-momentum tensor for the electromagnetic field in a Casimir cavity 
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placed in a weak gravitational field. The resulting calculations are considerably harder 
than in the case of scalar fields. By using Green-function techniques, we have evaluated 
the influence of the gravity acceleration on the regularized energy-momentum tensor of the 
quantized electromagnetic field between two plane-parallel ideal metallic plates, at rest in 
the gravitational field of the earth, and lying in a horizontal plane. In particular, we have 
obtained a detailed derivation of the theoretical prediction according to which a Casimir 
device in a weak gravitational field will experience a tiny push in the upwards direction [6] . 
This result is consistent with the picture that the negative Casimir energy in a gravitational 
field will behave like a negative mass. Furthermore, we find a trace anomaly in Eq. (4.7) 
proportional to the gravitational acceleration and vanishing for infinite plates' separation, 
not previously worked out for a Casimir device in a gravitational field. Our original results 
are relevant both for quantum field theory in curved space-time, and for the theoretical 
investigation of vacuum energy effects (see below). 

We stress that in our computation we do not add by hand a mass term for photons, 
unlike the work in Ref. [17], since this regularization procedure breaks gauge invariance 
even prior to adding a gauge-fixing term, and is therefore neither fundamental nor de- 
sirable [12,16]. In agreement with the findings of Deutsch and Candelas for conformally 
invariant fields [18], we find that on approaching either wall, the energy density of the 
electromagnetic field diverges as the third inverse power of the distance from the wall. It 
is interesting to point out that, in the intermediate stages of the computation, quartic di- 
vergences appear in the contributions from the ghost and the gauge breaking terms, which 
however cancel each other exactly. The occurrence of these higher divergences in such 
terms is also consistent with the results of Deutsch and Candelas, in view of the obvious 
fact that ghost fields are not ruled by conformally invariant operators. Unfortunately, a 
quantitative comparison with their results is not possible because they assume a traceless 
tensor, which is not the case in our problem where a trace anomaly is found to arise. 

Our results, jointly with the work in Refs. [6,19], are part of a research program aimed 
at studying the Casimir effect in a weak gravitational field, with possible corrections (albeit 
small) to the attractive force on the plates resulting from spacetime curvature [20] (cf. 
the recent theoretical analysis of quantum vacuum engineering propulsion in Ref. [21]). 
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Hopefully, these efforts represent a first step towards an experimental verification of the 
validity of the Equivalence Principle for virtual photons. 
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